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ABSTRACT 

EUV (Extreme-Ultraviolet) waves are globally propagating disturbances that have been observed since 
the era of the SoHO/EIT instrument. Although the kinematics of the wave front and secondary wave 
components have been widely studied, there is not much known about the generation and plasma 
properties of the wave. In this paper we discuss the effect of an EUV wave on the local plasma as 
it passes through the corona. We studied the EUV wave, generated during the 2011 February 15 
X-class flare/CME event, using Differential Emission Measure diagnostics. We analyzed regions on 
the path of the EUV wave and investigated the local density and temperature changes. From our 
study we have quantitatively confirmed previous results that during wave passage the plasma visible 
in the Atmospheric Imaging Assembly (AIA) 171 A channel is getting heated to higher temperatures 
corresponding to AIA 193 A and 211 A channels. We have calculated an increase of 6 - 9% in density 
and 5 - 6% in temperature during the passage of the EUV wave. We have compared the variation in 
temperature with the adiabatic relationship and have quantitatively demonstrated the phenomenon 
of heating due to adiabatic compression at the wave front. However, the cooling phase does not follow 
adiabatic relaxation but shows slow decay indicating slow energy release being triggered by the wave 
passage. We have also identified that heating is taking place at the front of the wave pulse rather than 
at the rear. Our results provide support for the case that the event under study here is a compressive 
fast-mode wave or a shock. 

Keywords: Sun: corona. Sun: evolution, waves 


1. INTRODUGTION 

Globally propagating disturbances in the solar atmo- 
sp here have been known sinc e they were first detected 
by |Moreton & Ramsey| (|1960|) in Ra images during flare 
observations, the so-called Moreto n wave s. Such phe¬ 
nomena were modeled by Uchidaj (1968) as the chro¬ 
mospheric imprint of fast mode IVlagnetohydrodynamic 
(MHD) waves in the corona. Though, it was not until 
the la unch of the Extreme-ultrav iolet Imaging Telescope 
(EIT, Delaboudiniere et al.|1995 ) c in-board the Solar and 
Heliospheric Observatory (SoHO, |Domingo et al.||1995) 
that the proposed coronal counterparts w ere identifled 
(jMoses et al. |1997t |Thompson et al.||1998|) . These tran- 
sient events have since been termed as “EIT waves” af¬ 
ter the instrument used to discover them. However, due 
to the debate regarding the true nature of these tran¬ 
sients and to incorporate observations made from differ¬ 
ent instruments, several authors refer to these events by 
other names s uch as “coronal bright fronts” (Gallagher 


^Lo^_2011), “coronal propaga ting fronts” ([bc^rijy^ 


et al.pufT fr^global EUV waves” (Patsourakos V Vourl- 
i das||2Ul2lh and “la rge-scale coronal propagating fronts” 
( [jNitta et al. 2013). In this paper we will refer to them 


kamalam.vanninathan@uni-graz.at 

^ Kanzelhohe Observatory/Institute of Physics, University of 
Graz, 9521 Treffen, Austria 


with a generic name such as EUV (Extreme-Ultraviolet) 
waves. 

Inco nsistencies in speeds between Moreton and EUV 
waves ( Klassen et al.|[2QQQ|, as well as the ob servation of 
stationary fronts (|Delannee & Aulanier|1999|) led to alter¬ 
nate interpretations and models for the observed wave¬ 
like signatures. There are two main opposing theories, 
waves and non-waves, to explain the visible features of 
EUV waves. In the non-wave models, these transient 
phenomena are interpreted as ground tracks of successive 
restructuring of magnetic held lin es during the eruption 


of a coronal mass ejection (CME, Delannee & Aulanier 
1999 Chen et al. 2002 Attrill et ar' ^0lj7| ). The wave 


plitude waves ( 

Thompson et al.||1998 Mann et al. 1999 

Warmuth et a 

[. 200ip. A third hybrid model tries to 


fronts: one consistent with the wave model and the o ther 
with the non-wave model (Zhukov & Auchere||2004). In 
recent year s, this theory h as gained support thro ugh ob¬ 
servations (Lin et al.|2010l) and MHD simulations (Cohen 


et al. 1|2009| iDowns et al.||2Ul‘2| 
TEe average speed ot EU V wc 


waves was measured to be in 


the range of 20 0 - 400 km s ^ in the EIT data (iThompson 
& Myers|2009|) and the Solar Terrestrial Relations Obser 


vatory (S'l'EREO) data ( Muhr et al.||2014 ) and updated 
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to 600 km s ^ from the Atmosp heric Imaging Assembly 


(AIA) data (Nitta et al. 2013). Fast EUV waves typi¬ 
cally decelerate during propagation, which has been in- 


fast magnetosonic waves ( 

jVrsnak et al.j 2006 

1 Lon^ 

y et al. 

20081 IVeronig et al.|2008 | 

Warmuth V Ivlann 

pHTT] 

IVluhr 

et al. 

|2014 

L Taking advantage of the quadrature conhgu- 


to study the 3D structure and evolution of EUV waves 


Pat sour akos et al. 


009[ Veronig et al 


H Kienreic^e^^ ]200^ 

[Ibrnmer et akpOTl 


Ma et al. 


2011). It was 

estiniated that most ot the EUV wave eriiission is com- 
i ng from heights of 80 - 100 Mm above the pho tosphere 


(jPatsourakos et al.||2009| |Kienreich et al.||2009|) which is 
comparable to 1 - 2 coronal scale heights (distance from 
the photosphere over which the pressure scales as a factor 
of 1/e). 

The high cadence and resolution of AIA has helped 
in improving the knowledge of kinematics of wave fronts 
and to identify secondary wave components. With re¬ 
gard to solar atmospheric research and plasma physics 
in general, it is hugely important to study how these 
waves affect the local plasma and pre-existing magnetic 
structures as they prop agate through. Oscillations and 
of filament s JjOkamoto et al.[2QQ4[ 

)nal loops (IWills-Davey & I'homi 


deflec tions _ 

2014), coronal loops ([Wills-Dav° 


Shen et al. 


. _ ^ Tiomp's Qn||l999|), 

coronal cavities (Liu et al.|2012[ ), and streamers ( jlripathi 


& Raouafi 


2007) have often been attributed to the pas¬ 
sage of eUV wave s. It has also been suggested that jets 
(iShen et al. || 2014 ) and sympathetic eruptions (Schrijver 


et al. can be triggered by such waves. 


4rom observations made with the EIT, the EUV Im¬ 
ager (EUVI) on-board the twin spacecrafts STEREO, 
and the AIA, EUV waves are typically seen as inten¬ 
sity enhancements in the 193 A (AIA), the 195 A (EIT, 
EUVI), the 211 A (AIA), the 284 A (EUVI), and some¬ 
times in the 335 A (AIA) channels while the 171 A chan¬ 
nel (EIT, EUVI, AIA) often shows co-spatial emission de¬ 
crease. This has prompted many authors to suggest that 
plasma from the 171 A channel is being heated, leading 
to the ionization of Ee ix/x to higher ionization states 
and which results in th e increased intensity as registered 
in the hotter channels (IWills-Davey & Thompsc 


Dai et al. 

2010 

bchrijver et al. |201ip 

range of t. 

he channels in which E U V ^ 


4'he temperature 

'waves are observed 

is about 1 - 2.5 MK. 

Studies on the plasma properties of EUV waves have 
been limited by the low chance of encountering such tran¬ 
sient phenomen a in the limited field -of-view (EOV) of 
spectrometers. Veronig et al.| ( 2Q11[ ) conducted a ded¬ 


icated observing campaign for registering EUV waves 
with the Hinode/EIS (EUV Imaging Spectrometer) in 
a sit-and-stare mode. The authors were successful in ob¬ 
taining one example where an EUV wave passed along 
the EOV of the EIS slit. The results from this campaign 
showed down-flows of 20kms“^ followed by up-flows of 



ferential Emission Measure (DEM) techniques tha t can 
be applied to imaging data. Kozarev et al. (2011) com¬ 


pared the wave front with a pre-event region using DEM 
analysis and noted a density increase of 12 - 18% at the 
location of the wave front. Their DEMs at the wave front 
also increased at higher temperatures which is consistent 
with plasma heating and compression. 

In this paper, we use DEM based analysis to study the 
response of the surrounding corona to the passage of a 
strong EUV wave associated with the first X-class flare 
of cycle 24 on 2011 Eebruary 15. The paper is organized 


as follows, in section [2 
of analysis. The resu 
discussed in section |4 


we describe our data and methods 
ts are presented in section and 
Einally, in section we give the 


summary and conclusions of this work. 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. Event overview 


The AIA (Lemen et al.|2012) instr ument on-board the 
Solar Dynarnics Observatory (SDO, Pesnell et al. |2012) 
images the Sun and its atmosphere using six EUV hl- 
ters. These filters put together cover a broad temper¬ 
ature range which facilitates the use of these data for 
the purpose of DEM analysis. In this paper we study 
the EUV wave associated with one of the strongest flare 
events that occurred after the launch of SDO/AIA. The 
2011 Eebruary 15 event is an X2.2 class flare within AR 
11158 (NOAA coordinates: S21 W28). This event has 
been the interest of previous studies relat ed to sun quakes 
( Kosovich e^ [2QjT Zharkov et al.||2Qll|), s unspot rota- 


tion (Jiang et al. 


2012[ Wang et al. '2014), cham »:es in 


the orient ation of phot ospheric magnetic fields (|Wang 
Gosain|2012), CME kinematics (iG^palswamy 


et al. 


et al. 


2012 


2012), (JME-CME int eraction ( TemmCT et al.|2014[ 


Shanmugaraju et al. 2014) and, o f particular relevance 


to the present stud y, EU V waves ( Schrijver et al.pOll 
Olmedo et al.||20^ b _ 


Lor this event, 


Schrijver et al. (2011) performed global 


magnetic field modeling and suggested that the observed 
EUV wave-like signatures are a result of compression 
between expanding coronal loops, which eventually be¬ 
come aCMEj_anRanoyer lying helmet streamer struc¬ 
ture. Olmedo et al. (2012) studied the interaction of the 
EUV wave with a coronal hole and detected reflected and 
transmitted wave components that are consistent with a 
fast mode MHD wave interpretation. 

This paper is related to the EUV wave that was asso¬ 
ciated with this event. The initiation time of the flare, 
according to GOES data, was 01:44 UT and the peak 
time was 01:56 UT. By examining the AIA 211 A channel 
images, we found that the first instance when the EUV 
wave pulse could be distinguished was at 01:50 UT and 
it was seen propagating isotropically across the disk. We 
studied the evolution of the EUV wave along two sectors 
which are at 55° and 210° (indicated by the red dashed 
lines in Eigure as measured in the anti-clockwise di¬ 
rection from the solar west. These particular directions 
were chosen because it was possible to observe a strong 
signal of the wave front, in these directions, for long dis¬ 
tances unhindered by other coronal structures such as 
active regions, coronal holes, coronal loops, etc. 

2.2. Perturbation profiles 
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Figure 1. First three columns show base-difference images of the evolution of the wave in 3 EUV channels (211 A, 193 A, 171 A) of AlA 
in which the wave is clearly visible. Fourth column shows the original data from 211 A channel. Co-temporal images are shown in each 
row. Time is marked at the bottom of each panel. The two sectors along which the EUV wave was studied are marked as 55 and 210, 
corresponding to the angles they make with the solar west as measured in the counter-clockwise direction, are shown enclosed by the red 
dashed lines. The green asterisks along each path marks the fixed ROI which was used to further investigate the plasma properties. 


As EUV waves cannot be easily visualized directly 
from original intensity images, images enhanced by a 
base-difference technique were used to study the abso¬ 
lute changes in the emission. The base image was taken 
at 01:45 UT and was subtracted from each subsequent 
image taken with one minute cadence. Although AIA 
data have a higher cadence of 12 s in the EUV channels, 
for the current study one minute cadence was sufficient. 
We only chose those images where the exposure times 
were constant and not affected by the automatic expo¬ 
sure control algorithm of the instrument. 

To be able to conveniently track the EUV wave as it 
propagates away from the CME site we made use of per¬ 
turbation profiles which are flux vs. distance plots to 
analyze the flux along the propagating path of the EUV 
wave. For calculating t he perturbation pro files we used 
the method adopted by Muhr et al. (2011). The helio¬ 
spheric coordinates of the flare given in the GOES data 
were assumed to be the center of the propagating wave. 
From this position we selected sectors of 5° angular width 
in the chosen directions. The propagation of the EUV 
wave is affected by local coronal structures and varies in 
different directions on the Sun. This study is related to 


the local changes in the quiet corona as the EUV wave 
passes through. Hence, we thought it necessary to choose 
a small sector so as to not average out all the properties 
of the EUV wave front. Along the chosen sectors, we 
defined successive annuli of 1° radial width between two 
constantly growing concentric circles from the wave cen¬ 
ter to study the properties of the EUV wave. The mean 
of all the pixel values in a given annulus within a sec¬ 
tor were plotted as a function of distance from the wave 
center to the limb. This procedure was repeated for each 
time step until the EUV wave faded away. The pertur¬ 
bation profiles were constructed from original images as 
well as the base-difference images. 


2.3. DEM analysis 

A DEM distribution gives information about the 
plasma distribution as a function of temperature along 
a given line-of-sight (LOS). For optically thin emission 
that is in thermal equilibrium a DEM, 0(T), is defined 
as 




(1) 
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where rie is the electron density at position h along the 
LOS at temperature T. For spectroscopic observations 
the line intensity, /a, is associated with the DEM by the 
equation 

h = J A{X) Gx{ne, T) 4>{T) dT, (2) 


2.4. Mean plasma temperature and density 


For this study, we made use of the HK code and de¬ 
rived the DEMs from the original da^ to calculate the 
DEM weighted average temperature (T) and plas ma den¬ 


sity (n) usi ng the following expressions (defined in Cheng 


where A{X) is the elemental abundance with respect to 
hydrogen and Gx{ne^ T) is the contribution function for a 
given spectral line. Analogous to this, for broad/narrow 
band filter observations the intensity in a filter, Ifuter^ is 
associated to the DEM by the equation 


^filter — j RfilteriT) 0(T)dT, 


(3) 


where RfuteriT) is the temperature response function of 
the instrument. With the help of inversion techniques 
or forward fitting models, the DEM from multi-channel 
o bservations can be derived using this relationship. 


Hannah & Kontar (2012) have developed a regularized 
inversion method to reconstruct the DEM from the six 
AIA channels that are sensitive to coronal temperatures. 
This method has been tested for various cases and was 
found to be reasonably robust and comp utationally fast 


(Del Zanna 2013 

their code provides error bars in 


[Plowman et al. 


2013|) . Additionally, 
3oth vertical (DEM) 


and horizontal (temperature) directi ons. However, this 
metho d is also with its limitations. IHannah fc Kontarl 
(|2013|) have cautioned about the DEM solution at regions 
with large uncertainties (>0.3 dex) in temperature res¬ 
olution. The application of this method also shows some 
spurious high temperature emission which is most likely 
contributions of cool lines to the AIA 94 A and 131 A 
channel s that are currently unacc ounted for in the atomic 
models ( Young fc Muglach|2014 ). The phenomenon that 
we study in this paper is known to be formed within the 
thermal range of about 1 MK to 2.5 MK. The above men¬ 
tioned caveats are not affecting this temperature range. 
Hence for the current study we have applied this DEM 
technique (henceforth referred to as HK code). 

Average counts in each of the six EUV filters of AIA 
were determined from the perturbation profiles and used 
as an input for the DEM code. Apart from data num¬ 
bers (DNs), the DEM code also accepts errors in DNs 
as an input. To cal culate the uncertainties, we used the 
method outlined in Yuan fc Nakariakov (|2Q12) for AIA 
data, which has been modified from |Ascnwanden et al 


(2000) for the use with TRACE data. The square of the 
total noise {crnoise) was determined using the formula 

^noisei^) ^photoni^) Y ^readout Y ^digit 

^compress ^dark ^subtract 


a 


spikes 


(F) 


which corresponds to the sum of the squares of un¬ 
certainties, in units of DNs, in photon Poisson noise 
{cr photon), electronic readout noise {cr readout), digitiza¬ 
tion noise {cr digit), compression noise {cr compress), dark 
current noise {(Jdark), subtraction noise {cr subtraet) , noise 
due to removal of spikes in the images {crspikes), and er¬ 
rors in the response functions {(Jresv )- The values for the 
i nstru ment al errors were taken f rom I Yuan fc NakariakovI 
(2012) and Boerner et al. ( 2Q12| ). 


T = 


n = 


f c/)(T)TdT 
imdT ’ 

(5) 

^IlmdT 

(6) 


where h is the distance along the LOS. From quadrature 
studies, the height of the EUV wave s have been estab¬ 
lished to be between 80 - 1 00 Mm (Patsourakos et al. 
2009 Kienreich et al. 2009). For the calculations here 
we assume h = 90 IVlm. 

We also made difference DEMs (referred to as ADEM 
henceforth) by subtracting the DEM of the base image 
(taken at 01:45 UT) from the DEM at each subsequent 
time step. This was used to better illustrate the changes 
in the DEM during the passage of the EUV wave. This 
can be considered to be similar to background subtracted 
images or base images for the purpose of image enhance¬ 
ment. 


3. RESULTS 


800 



Ol_^^^^^^^^^^^^^^^^^^^_ 

01:56 02:00 02:04 02:08 02:12 

Start time (15-Feb-11 01:53:00 UT) 

Figure 2. Distance vs. time plot along with Icr errors obtained 
from AIA 211 A channel for the propagation of the EUV wave front 
along 55° and 210° direction. 


In order to study the response of the corona to the 
passage of the EUV wave along our chosen directions of 
study, we selected regions of interest (ROIs), hereafter re¬ 
ferred to as ROI55 and ROI210 corresponding to 55° and 
210° from the solar west, respectively, (marked by green 
asterisks in Figure and studied the plasma properties 
of these regions as a function of time. The size of ROI55 
was chosen such that at a given time step (01:58:48 UT 
for this direction) the whole wave pulse occupies the ROI. 
This size was found to be ~ 168 Mm. The same could not 
be done for ROI210 due to the presence of a stationary 


brightening which is discussed later in Section 3.2 


Figure ^ shows the EUV wave kinematics in the 
AIA 211 A channel along the two sectors covering the 
ROIs. The position of the EUV wave (indicated by the 
blue arrows in Figures]^ andwas determined from the 
perturbation profiles of the 211 A channel, as the wave is 
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R0I55 



Figure 3. Perturbation profiles along the 55° direction derived 
from AIA 211 A base-difference images at five different time steps 
as indicated in the panels. The blue arrow points to the prop¬ 
agating EUV wave pulse. The section of the plots indicated in 
red corresponds to R0155 marked by green asterisks along the 55° 
direction in Figure 

clearly distinguishable here, and the same position was 
used in the remaining channels. We identified the leading 
front by fitting a Gaussian function to the wave pulse. 
The MPFIT routine was used to calculate the Gaussian 
parameters along with its errors. The leading front, xiead, 
was defined as 

a;iead =Pl +2p2, (7) 

where pi is the centroid of the Gaussian and p 2 is the 
width. Accordingly, the la errors in xiead is given by 

Aa^iead = ± (Api + 2Ap2). (8) 

The speed of the propagating wave front thus estimated, 
in the 211 A channel, along the 55° direction was found 
to be 437± 19kms“^ and along the 210° direction was 
366± 44kms“^. 

3.1. Study along 55° direction 

In Figurewe show the perturbation profiles from AIA 
211 A channel along 55° direction at five different time 
steps corresponding to before (panels 1,2), during (panel 
3) and after (panels 4, 5) the passage of the FUV wave 
from ROI55, marked in red. The average flux from this 
region taken from all the six AIA coronal channels were 
used as inputs for the DFM analysis. 

The curves resulting from our DFM analysis are shown 
in Figure The left panels correspond to DFM curves 


derived from the perturbation profiles constructed from 
the original data. The right panels are ADFM curves de¬ 
rived by subtracting the DFM of the base image (taken 
at 01:45 UT) from the DFMs shown in the left panels. 
Here, the negative values imply a reduction in emission 
with respect to the pre-event corona. These curves are 
plotted for the same time steps as the perturbation pro¬ 
files shown in FigureThe DFMs in the left panels were 
used for further calculations while the ADFMs in the 
right panels illustrate the temperatures at which there 
is a reduction/increase in emission compared to the pre¬ 
event corona. These changes are not very obvious from 
the original data DFMs (left panels). 

In Figure the top two rows show the DFM and 
ADFM before the FUV wave arrived at ROI55. The 
right panels show very small change in DFM values and 
is comparable to the background emission. The middle 
row corresponds to the DFM and ADFM when the FUV 
wave pulse was at ROI55 while the bottom two rows of 
this figure shows the DFM and ADFM when most/all of 
the wave pulse has moved out of ROI55. 

Using equations and we derived the averaged tem¬ 
peratures and densities, respectively, from the DFMs ob¬ 
tained with the original data (left panels of Figure]^. 
Figure shows the time evolution of the derived quanti¬ 
ties along with its errors. The temperatures and densities 
obtained in the time range of 01:50 UT and 01:53 UT are 
approximately constant as the FUV wave had not yet 
reached ROI55. We took the average of these quanti¬ 
ties as the background temperature and density. From 
the values plotted in Figure we measured an 8.7% 
increase in density compared to the background (from 
1.87xl0^cm“^ to 2.03x 10^ cm“^) as the FUV wave 
passed through ROI55. A similar calculation for temper¬ 
ature showed a 6.3% increase from 1.48 MK to 1.57 MK. 

After the FUV wave left ROI55 we saw that at 
first the density falls back to the background levels (at 
02:06:48 UT in Figure and then it increases slightly 
until 02:12:00 UT. However, during this time the tem¬ 
perature is not seen to reach the background levels and 
it remains at an elevated level. 

We calculated the absolute errors in density to be of 
the order of 0.16xl0^cm“^ and those in logioT to be 
of the order of 0.06 dex. We note here that the main 
contribution to the error in density is from the error in 
the value of h which is ~10%. In this paper we discuss 
about the relative changes and not the absolute values. 
Hence, we want to determine errors which are more rele¬ 
vant to the current study rather than using the absolute 
errors. We have shown that the temperature of the event 
studied here varies between logicT=6.17 and 6.19. This 
falls within the range where the AIA channels have very 
good temperature sensitivity. So it is more appropriate 
to discuss errors for changes within this specific tempera¬ 
ture range rather than for the whole. Synthetic Gaussian 
DEM curves were folded through the AIA instrument re¬ 
sponse and the result was used as input to the HK code 
to obtain DEM reconstructions. The mean temperature 
and total emission measure calculated using these curves 
were used to derive the errors. A detailed description 
of the method is given in the appendix. The errors ob¬ 
tained, using this method, in density is of the order of 
0.014 xlO^cm ^ and in temperature is 0.004 dex. Fig- 
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Originol DEM 



Figure 4. DEM curves from the original images (left) and ADEM 
data DEMs (right) for ROI55 at the time steps indicated in each o 


Base-difference DEM 



obtained by subtracting the base DEM (at 01:45 UT) from the original 
f the panels. Same as the time steps given in Figure]^ 
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Figure 5. Average density (top) and temperature (bottom) de¬ 
rived from the DEM evolution at ROI55 as a function of time. Thin 
black lines represent the absolute errors while the thick shorter lines 
show the relative error. 


ure[^ shows for comparison bo th, t he a bsol ute errors and 


errors derived using equations |A3| and A4[ In the subse¬ 
quent images we show only the newly derived errors that 
relate to the relative changes in density and temperature. 
To verify the theory of adiabatic compression of the 
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Figure 6. The expected change in temperature for an adiabatic 
process are plotted as blue plus symbols and the temperature 
changes obtained from the DEM diagnostics are plotted as red 
diamonds. The absolute temperature scale is plotted on the right 
axis. For ease of viewing the errors are shown in the legend rather 
than on the data points. 

plasma at the wave front we used the adiabatic relation 
between temperature and density 



(assuming the adiabatic index 7 = 5/3 for fully ion¬ 
ized plasma). We calculated the change in temperature 
(T/Tq) from our DEM estimated values and compared 
them with the expected change in temperature during 
an adiabatic process as derived from the density changes 
(p/po) using equation]^ (see Figure]^. We found that 
during the heating phase the observea and expected val- 
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ues are in good agreement but the same is not true for 
the cooling phase. We also note that the secondary slow 
increase in density from 02:06:48 UT to 02:12:00 UT is 
not resulting in an adiabatic increase of temperature. 

We did the same analysis as described above by de¬ 
creasing the size of ROI55 to 96 Mm. In this case of finer 
spatial binning we see that the general evolution is the 
same but the derived changes are larger: a 9.9% increase 
in density and a 6.1% temperature increase with respect 
to the background values (see Figure [^. Again, the heat¬ 
ing phase is consistent with adiabatic compression at the 
wave front while the cooling phase is not. 


small R0I55 



01:52 01:56 02:00 02:04 02:08 02:12 

Stort Time (15-Feb-ll 01:49:00 UT) 


Figure 7. Sample plots for test done with finer spatial binning at 
ROI55. Average density (top) and temperature (middle) derived 
from the DEM evolution as a function of time. Bottom: The 
expected change in temperature for an adiabatic process are plotted 
as blue plus symbols and the temperature changes obtained from 
the DEM diagnostics are plotted as red diamonds. The absolute 
temperature scale is plotted on the right axis. For ease of viewing 
the errors are shown in the legend rather than on the data points. 


3.2. Study along 21 (F direction 

For comparison we performed a similar study along the 
210° direction. We show example plots of perturbation 
profiles in Figure with the position of the FUV wave 
indicated by the blue arrows. In this direction the clear 
identification of the FUV wave was complicated by the 
presence of a stationary brightening close to the location 
of the flare. The position of the brightening is indicated 
by a horizontal green dashed line in the panels of Figure 
where it appears. Such bright enings have been previously 
eported for different events (iDelannee fc Aulanier||1999 


Attrill et al.||2QQ7 Cohen et al.||2UU9| ~ lVluhr et al.||2Ull 


From Figure we see that the brightening does not 
exist prior to the passage of the FUV wave (top panel). 
By comparing panel 2 with panels 3, 4 and 5 we can say 
that the brightening occupies the same position as the 
first appearance of the FUV wave pulse and continues to 
remain there at least for the next 15 min. In panel 5 of 
this figure we see that a secondary stationary brightening 
is formed adjacent to the first one after the FUV wave 
pulse has moved forward. The clear distinction of the 
wave pulse from the stationary brightening is difficult. 
We selected ROI210 such that it does not include the 
primary stationary brightening. 

W e fo llowed the same procedure as described in Sec¬ 
tion 3.1 and constructed DFM curves from the original 
data for ROI210. We used the DFMs to determine den¬ 
sities and temperatures and studied the time evolution 
of these quantities (see Figure]^. 
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Figure 8. Perturbation profiles along the 210° direction derived 
from AIA 211 A base-difference images at five different time steps 
as indicated in the panels. The blue arrow points to the propa¬ 
gating EUV wave. The section of the plots indicated in red cor¬ 
responds to RO1210 marked by green asterisks along the 210° di¬ 
rection in Figure^ The dashed green horizontal line marks the 
position of the stationary brightening. 


From Figure we found that the relation between den¬ 
sity and temperature is different from what we observe 
in ROI55. The density in ROI210 peaks at 02:02 UT, the 
temperature continues to rise for 4 mins after this time 
until 02:06 UT before it begins to fall. From the val¬ 
ues plotted in Figure we calculated a 5.7% increase in 
density and a 4.5% increase in temperature. The density 
increased from 1.97x10^ cm“^ to 2.08x10^ cm“^ and the 
temperature increased from 1.36 MK to 1.42 MK. 
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R0I210 
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Figure 9. Same as Figure but for ROI210. 

We also calculated the adiabatic changes for ROI210 
and found it to be similar to ROI55 with one exemption. 
From Figure!^ (bottom) we see that the temperature rise 
phase until 02:02 UT nicely fit with the adiabatic rela¬ 
tionship, but the observed increase thereafter, 02:02 UT 
to 02:06 UT, is not consistent with heating due to adia¬ 
batic compression and requires other additional agents. 

4. DISCUSSIONS 

The appearance of the EUV wave as a dark front in 
the 171A channel has often been interpreted as emission 
decrease in this channel due to heating. In our study 
with the ADEM curves we have shown that, when the 
EUV wave passes over the ROI there is an increase in the 
DEM values in the region bounded by 6.1 < logio T < 6.4 
(which coincides with the peak of the temperature re¬ 
sponse curves for the AIA 193 A and 211A channels) 
while there is a decrease in the DEM in the region of 
5.8 < logio T < b.l (corresponding to the peak temper¬ 
ature of the AIA 171 A filter). Erom our observations 
the average pre-event temperature Tq = 1.48 MK (see 
Eigure and lies to the right of the response peak of 
the 171 A channel (along the negative gradient of the re¬ 
sponse curve) and to the left of the response peaks of the 
193 A and 211 A channels (along the positive gradient of 
the response curve). This implies that the mean DEM- 
weighted temperature is shifting locally along the re¬ 
sponse curve, R(T), and the change in emission is propor¬ 
tional to the derivative of the response curve dR(T)/dT. 
A temperature increase from Tq means it is moving down 


along the 171 A response curve and up along the 193 A 
and 211 A response curves, implying a decrease in flux in 
the 171 A channel (darkening) and an increase in flux in 
the 193 A and 211 A channels (brightening). This con¬ 
firms that the dark front observed in 171 A channel is 
a result of plasma being heated to higher temperatures. 
After the EUV wave has passed the ROI, we notice that 
the peaks of the ADEM values are diminished and have 
shifted slightly towards lower temperatures indicating 
that there is cooling of the plasma behind the wave. 

In the relaxation phase at ROI55, we have pointed out 
that the density falls back to background levels. In case 
of a coronal dimming behind the wave we would expect 
the density to drop below the background values. This 
observation suggests that the EUV wave is not followed 
by a dimming region or in other words, this is a wave phe¬ 
nomenon which only displaces the plasma in its current 
position and does not carry it forward. 


AIA 211 channel 02:04:48 UT 
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Figure 10. Top: Perturbation profile from the 2IlA channel 
taken at 02:04:48 UT. The vertical lines divide the wave pulse into 
three segments: front, peak and rear (from right to left). Bottom: 
The DEMs for the three segments shown above and as indicated 
by the legend. 

We have also studied in detail the plasma characteris¬ 
tics across the wave pulse by dividing it into three seg¬ 
ments corresponding to front, peak and rear as shown in 
Eigure [Tq| We found that for a particular time step the 
segments corresponding to the peak and front of the pulse 
were hotter and denser than the rear. This sequence was 
true for all the time steps where a clear identification of 
the wave pulse was possible. Average temperatures for 
the front, peak and rear segments of the wave pulse were 
1.75T0.008MK, 1.77T0.008MK and 1.60T0.007MK, re¬ 
spectively. Such an observation indicates that heating 
due to compression is taking place at the wave front 
which is consistent with the behavior of a compressive 
wave. 

The relaxation phase in ROI55 shows that the tem¬ 
perature remains at an elevated level while the density 
returns to the background levels. This deviation could be 
due to non-linear effects at play for a such a global wave 
pulse, which are not accounted for in equation The 
peak density amplitude derived is 6%, which is actually 























Plasma diagnostics of an EUV wave 


9 


quite big considering the assumption of small perturba¬ 
tions (linear theory) from which the standard magneto¬ 
acoustic waves are derived. Alternatively, it could be 
a result of local energy release triggered by the pass¬ 
ing wave fr ont, e.g. reconnectio n with favorably oriented 
field lines (Attrill et al. 2007), which provides heating 
in addition to the adiabatic process. There have been 
several examples (as mentioned in Section of EUV 
waves disturbing structures that are present along its 
path. When the EUV wave passes over the quiet-Sun, it 
disturbs the local magnetic field lines which can result in 
magnetic reconnection behind the wave front. Although 
reconnection is not the primary mechanism producing 
the wave front, it can be a side effect of the passing wave. 
Such occurrences will contribute in increasing the local 
temperature while it need not necessarily change the lo¬ 
cal density. 

Based on the relative emission changes at the wave 


front in the A IA 171 A, 193 A and 211 A filters, Schrijver 


et al. (2011) estimated that the observations are con¬ 
sistent with adiabatic compression for plasma tempera¬ 
tures between 1.2 and 1.8 MK. Assuming adiabatic com¬ 
pression as the only source of heating, they have shown 
that the observed intensities in the three filters are con¬ 
sistent with adiabatic compression at the wave front for 
T=1.3MK. Assuming isothermality, they derived a max¬ 
imum density increase of 10% associated with a temper¬ 
ature increase of 7% for the east-west propagation of the 
EUV wave over the quiet-Sun. They interpreted these 
enhancements as a result of compression of the expanding 
CME against the stalks of the overlying helmet-streamer. 

We have taken a different approach and have per¬ 
formed a DEM analysis to study the plasma properties 
at the wave front. Our study is a detailed quantita¬ 
tive analysis of the compression and relaxation phase, 
at and behind the wave front, respectively. The in¬ 
crease in density and temperature of 9% and 6%, respec¬ 
tively, in the 55°direction and 6% and 5%, respectively in 
the 210°direction that we have calculated in the current 
study are in basic agreement with the results of Schrijver 
et al. (2011). However, we interpret this phenomenon as 
a compressive fast-mode wave for the following reasons: 
(1) we have found that the evolution of temperature and 
density are consistent with adiabatic compression, as is 
expected in the case of a compressive wave; (2) we are 
able to detect a wave pulse quite close to the active re¬ 
gion while in the streamer scenario the increase in in¬ 
tensity should only be visible when the CME approaches 
the streamer footpoints; (3) we have analyzed different 
segments across the wave pulse and found that heating 
and compression is taking place close to the front of the 
wave which is the likely scenario for a compressiv e wave 


whereas in the c urrent shell model (as suggested in Schri- 
jver et al.||2011|), the highest temperature is expected to 
be close to the rear part of the pulse which is adjacent 
to the current layer. 

The stationary brightening that we detect in the 210° 
direction has been previously interpreted as compression 
of plasma near the foot points of opening field lines due 


to the expansion of CM E flanks (Delannee & Aulanier 


1999 Muhr et al. 2011) or as ongoing reconnection be¬ 


tween expanding flux r opes and surround i ng favorably 


oriented magnetic field (Attrill et al. 2007 Cohen et al. 


2009). In either case this could result in energy release 
and associated temperature increase. We see the effect 
of this when we compare the densities and temperatures 
at ROI210. Firstly, the density in ROI210 does not re¬ 
turn to the background level as it did in ROI55. Sec¬ 
ondly, there is an increase in temperature even after the 
density has reached its peak. These observations can 
be attributed to the presence of the secondary station¬ 
ary brightening which occurs at ROI210 (see panel 5 of 
Figure!^. We also confirm this from the comparison of 
actual temperature change with the expected adiabatic 
change. The temperature increase from 02:02 UT until 
02:06 UT is not due to adiabatic compression. 

In our study we show that the theory of adiabatic com¬ 
pression at the wave front is consistent for both directions 
that we have studied, with the exclusion of the effect of 
the stationary brightening in the 210° direction. How¬ 
ever, the relaxation phase is more difficult to account 
for. Since the EUV wave alters structures on its path, 
it cannot be expected that a particular region remains 
the same before and after the passage of the wave. De¬ 
viation from the adiabatic relationship in the relaxation 
phase could be attributed to local changes triggered by 
the passage of the EUV wave that provide additional re¬ 
lease of energy through a process which is not adiabatic. 


5. SUMMARY AND CONCLUSION 

We present here a study on the EUV wave associated 
with the large X-class flare of 2011 February 15. We 
studied the response of the quiet corona to the passage 
of the EUV wave using DEM analysis. We selected two 
different directions to study the wave. One was “undis¬ 
turbed” by external factors (at 55°) while the second was 
“disturbed” by the presence of a stationary brightening 
(at 210°). 

In the 55° direction, we have shown that the wave front 
approaching ROI55 compresses the plasma underneath 
resulting in heating that is purely adiabatic in nature. 
There is a 6% increase in temperature corresponding to 
a 9% increase in density. When the EUV wave departs 
from ROI55 we see that the density falls back to val¬ 
ues close to its initial state suggesting that there is no 
dimming region behind the EUV wave. This observation 
supports that the phenomenon is a wave only displac¬ 
ing plasma in its current position. On the other hand, 
the temperature drops slowly and does not return to its 
background values. It remains at an elevated level after 
some time. 

In the 210° direction, the approaching EUV wave front 
at first compresses the plasma adiabatically. 5% increase 
in temperature is associated with 6% increase in density. 
Later a stationary brightening is formed at the position 
of ROI210 which affects the density and temperature of 
the local plasma. We notice that when the EUV wave 
has departed ROI210 the temperature continues to rise 
briefly despite the density decrease indicating that the 
wave passage has induced some local energy release pro¬ 
cesses. Unlike ROI55, the density in this region does not 
return to its initial state. 

In this study, it has been quantitatively shown that 
the process of adiabatic compression is heating the coro¬ 
nal plasma at the wave front, consistent with a MHD 
fast mode wave. The behavior of the EUV wave dur¬ 
ing the heating phase, as we have shown here, provides 
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support to the wave nature of the phenomenon under 
study. The temperature structure across a wave pulse 
shows heating is taking place at the front of the wave 
and this provides additional evidence of the event being 
a fast magnetosonic wave. 
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APPENDIX 


We have used the following method to test the sensi¬ 
tivity of the HK code that is used in this study and to 
derive errors within the temperature range under consid¬ 
eration here. A typical DEM curve is most similar to a 
Gaussian curve. So we have used a Gaussian, drawn as 
a function of temperature, as a proxy for a DEM. The 
amplitude and width of the Gaussian is taken to be close 
to the observed DEMs. Since the temperature of the 
event studied in this paper varies from logicT=6.17 fo 
6.19 for the 55°direction and logioT=6.13 to 6.15 in the 
210°direction, we varied the position of the center of the 
Gaussian peak from logioT=6.10 to 6.30 in intervals of 
0.01. The temperature sensitivity of all the AIA EUV 
channels are the best within the temperature range of 
the event studied in this paper. So we use this method 
to derive errors that is most relevant to this study, where 
we are mostly interested in the changes in plasma tem¬ 
perature and density (rather than the accuracy of the 
absolute values). Using the relation 


I filter = ^Rfilter{T) f{T) AT, 


(Al) 


where RfuteriT) corresponds to the temperature re¬ 
sponse function of an EUV filter in AIA and /(T) is 
the input Gaussian representing a DEM, we were able 
to obtain the simulated counts for the six EUV chan¬ 
nels. These counts were further used as inputs for the 
HK code to obtain the corresponding DEMs. A compar¬ 
ison between the input Gaussian and the DEM from HK 
code is shown in Eigure IE Erom these DEM curves, 
the weighted average temperature and emission meas ure 
[EM(T)) were calculated using equationsl^ and A2, 


spectively and plotted as a function of the Gaussian peak 
temperature: 


EM{T) = / 


dT. 


(A2) 


This exercise was repeated by introducing 3 different lev¬ 
els of random errors (3%, 5% and 10%) to the input 
Gaussian curves that represent DEMs. The results from 
introducing 5% error is shown in Eigurej^ Eor each level 
of error applied to the input Gaussian, 100 realizations 
of randomness were performed. Erom the 100 realiza¬ 
tions, the mean and standard deviation of the emission 
measure and temperature were calculated for each error 
level. These plots for the 5% random error are shown 
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Figure 11. An example of an input Gaussian curve along with 
5% error (top) and its corresponding DEM curve (bottom). 
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Figure 12. The emission measure (top) and temperature (bot¬ 
tom) as calculated from the HK code. The dashed lines in each of 
the plots depicts the values calculated from the input Gaussian. 

in Eigure We notice that the HK code is able to 
reproduce the Gaussian input temperatures quite accu¬ 
rately with a standard deviation of 0.0004 dex for 5% 
error while the standard deviation for emission measure 
is 2.33x10^^ cm“^ in the considered range of Gaussian 
peak temperature (logioT=6.10 to 6.30). In this study 
we are interested in the changes in the plasma parameters 
so we want to derive errors which are more relevant than 
just the absolute errors. We use the following relations 
to define the typical errors in temperature and emission 
measure within the temperature range under study here: 

AEM = \EMdem — EMinput\, (A3) 


AlogT = \logTDEM — log 

Einput I, (A4) 

where T^eMi EMdem correspond to the values plot¬ 
ted in Eigure [T^ and Ti^put^ EMinput correspond to the 
dashed lines shown in the same figure. Using this method 
we obtain errors in temperature as 0.004 dex and in emis¬ 
sion measure as 5.2x10^^ which corresponds to an error 
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of 0.014 xlO^cm ^ in density. These were used instead 
of the absolute errors since they are more specific to 
the temperature range associated with the event we are 
studying. Comparison between the two errors is shown 
in Figure 
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